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ABSTRACT 

Context. The properties of the early-type binary Cyg OB2 #5 have been debated for many years and spectroscopic and photometric 
investigations yielded conflicting results. 

Aims. We have attempted to constrain the physical properties of the binary by collecting new optical and X-ray observations. 
Methods. The optical light curves obtained with narrow-band continuum and line-bearing filters are analysed and compared. Optical 
spectra are used to map the location of the He ii A 4686 and Ha line-emission regions in velocity space. New XMM-Newton as well 
as archive X-ray spectra are analysed to search for variability and constrain the properties of the hot plasma in this system. 
Results. We find that the orbital period of the system slowly changes though we are unable to discriminate between several possible 
explanations of this trend. The best fit solution of the continuum light curve reveals a contact configuration with the secondary star 
being significantly brighter and hotter on its leading side facing the primary. The mean temperature of the secondary star turns out 
to be only slightly lower than that of the primary, whilst the bolometric luminosity ratio is found to be 3.1. The solution of the light 
curve yields a distance of 925 ± 25 pc much lower than the usually assumed distance of the Cyg OB2 association. Whilst we confirm 
the existence of episodes of higher X-ray fluxes, the data reveal no phase-locked modulation with the 6.6 day period of the eclipsing 
binary nor any clear relation between the X-ray flux and the 6.7 yr radio cycle. 

Conclusions. The bright region of the secondary star is probably heated by energy transfer in a common envelope in this contact 
binary system as well as by the collision with the primary's wind. The existence of a common photosphere probably also explains 
the odd mass-luminosity relation of the stars in this system. Most of the X-ray, non-thermal radio, and possibly y-ray emission of 
Cyg OB 2 #5 is likely to arise from the interaction of the combined wind of the eclipsing binary with at least one additional star of this 
multiple system. 

Key words. Stars: early-type - Stars: individual: Cyg OB2 #5 - Stars: binaries: eclipsing - Stars: winds, outflows - Stars: fundamental 
parameters - X-rays: stars 



1. Introduction 

; The eclipsing nature of CygOB2 #5 (= V 729 Cyg = 
! BD -h40° 4220) was first reported by Miczaika (W53). 
Photometric data were subsequently published by Hall (1974), 
who gathered 143 UBV observations. Hall (1974) noted the 
existence of intrinsic variability as well as a slight asymmetry 
in the light curve which made the system appear brighter during 
the maximum that followed the secondary eclipse. The light 
curve obtained by Hall ( |1974| l revealed unequal eclipse depths. 
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thus suggesting that the primary should have a larger surface 
brightness than the secondary. 

From the analysis of a set of optical spectra, Bohannan & 
Conti d 19761 1 concluded that the two components of Cyg OB2 
#5 should be of equal brightness, whilst they simultaneously in- 
ferred a mass ratio (primary/secondary) of ^ = 4.3. Therefore, 
the secondary star would appear to be far too luminous for 
its mass, and these authors suggested accordingly that the sec- 
ondary might be an O-star on its way to becoming a Wolf-Rayet 
object. A re-analysis of the same spectra by Massey & Conti 
(119771 ) yielded a somewhat lower mass ratio of 3.3, but did not 
alter the basic conclusions of the Bohannan & Conti d 19761 1 pa- 
per 

Leung & Schneider ( fT978b analysed Hall's ( fT974l l UBV light 
curve adopting the mass ratio proposed by Bohannan & Conti 
( 119761 ). Their description of the system was that of an overcon- 
tact binary with the primary accounting for almost 90 % of the 
total light, unlike the almost equal luminosities inferred from the 
strengths of the He absorption lines (Bohannan & Conti [T976] l. 
Whilst the luminosity difference of 2. 1 mag inferred by Leung & 
Schneider ( 1978 ) is in qualitative agreement with a large mass 
ratio, it is clearly at odds with the strong spectral signature of the 
secondary in the spectrum of the binary. 
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Vreux (I1985I I studied the variations in the Ha emission Hne 
profile and suggested that the primary is currently losing mate- 
rial towards the secondary. This was interpreted as further evi- 
dence that the secondary star could be a Wolf-Rayet type object 
with an unusual spectrum produced by the accretion of the ma- 
terial. 

Rauw et al. (1999') analysed an extensive set of optical spec- 
tra and presented a revised orbital solution along with a discus- 
sion about the phase-locked spectroscopic variability. The new 
orbital solution yielded a mass ratio of ^ - 3.31 + 0.25 and 
revealed a large difference in the systemic velocities of the pri- 
mary and the secondary (the latter having a more negative y). 
Based on the relative strengths of classification absorption lines, 
Rauw et al. ( 119991 ) assigned an 06.5-7 spectral type to the pri- 
mary and Ofpe/WN9 to the secondary. They also compared the 
strengths of the primary and secondary absorption lines to infer 
a visual brightness ratio of 1.4 + 0.6. The complex behaviour of 
the absorption lines (some of which go into emission during part 
of the orbit) prevents however a more accurate determination of 
this ratio. The radial velocities of the emission lines imply that 
the emissions most likely originate in material that is moving 
from the primary towards the secondary. Rauw et al. d 19991 1 also 
reported on variations in the equivalent width of the He n 4686 
emission line, which suggest an occultation effect of the emis- 
sion zone in addition to the modulation of the continuum light 
curve due to the photometric eclipses. The variations in the line 
EW suggest that the emission arises from a region close to the 
secondary star, slightly preceding this star on its orbit. One of 
the main motivations of the present photometric campaign was 
to study this behaviour using a set of specially designed narrow- 
band photometric filters. 

Another goal of this paper is to study the variability of the 
X-ray emission of CygOB2 #5. X-ray emission from this star 
was first discovered with the EINSTEIN satellite in December 
1978 during a calibration observation of CygX-3 (Harnden et 
al. |1979l ). Since then, our target has been observed with many X- 
ray satelUtes, including /;C»5Ar (Waldron et al. [T998] l and A5CA 
(Kitamoto & Mukai 1996). We conducted an observing cam- 
paign with XMM-Newton to improve our knowledge of the X-ray 
emission of Cyg OB2 #5. 

Finally, we emphasize that important information about this 
enigmatic system can be derived from the radio domain. Indeed, 
Miralles et al. ( fT994l l noted that the radio emission of Cyg OB2 
#5 alternates between a low state with a spectral index that is 
roughly consistent with thermal wind emission and a high state 
with a much flatter (partially non-thermal) spectral index on a 
~ 7yr period. Abbott et al. (|1981| l and Mirafles et al. (1994) 
reported on the existence of a radio companion 0.8 arcsec to the 
NE of the main radio source (which is associated with the eclips- 
ing binary). Subsequent observations by Contreras et al. d 19971 1 
revealed that this radio companion has an elongated shape and 
lies in-between the short-period binary and a third stafl which 
was first reported by Herbig ( |1967l l. Contreras et al. d 19971 1 sug- 
gested accordingly that the radio companion corresponds to the 
wind interaction zone between the binary system and the ter- 
tiary component. Recently, Kennedy et al. (2009) re-analysed all 
VLA observations of Cyg 0B2 #5 and showed that the primary 
radio source, associated with the eclipsing binary varies on a pe- 
riod of 6.7 + 0.2 yr whilst the secondary source remains constant. 



' The latter star is 3 - 4 mag fainter than the binary system, and its 
colours and magnitude suggest a B0-B2V spectral type. The angular 
separation between this astrometric companion and the close binary 
system suggests an orbital period of several thousand years. 



The variations in the main radio component can be represented 
by a simple model in which a fourth (unresolved) object orbits 
the eclipsing binary in an eccentric orbit and the varying radio 
emission results from a variable non-thermal emission produced 
in the wind interaction between the short period binary and the 
fourth star (Kennedy et al. |2009l l. 

2. Observations 

2.1. Photometry 

Our photometric monitoring campaign was conducted over 14 
nights in August - September 1998 (between HJD 2451047.3 
and 2451060.6) at the OHP 1.20 m telescope. The telescope was 
equipped with a 1024 x 1024 SITe CCD camera. The pixel size 
was lAfim corresponding to 0.69 "on the sky. We used a set of 
five filters with narrow bandpasses designed especially for the 
photometric search and classification of Wolf-Rayet stars (see 
Royer et al. |1998l l. Three of these filters are centred on lines that 
are in emission in various subtypes of Wolf-Rayet objects (He ii 
.14686, Hei .15876, and Civ AA5?,Q\, 5812), whilst the other 
two (hereafter ci and C2 respectively) are centred on line-free re- 
gions around 5057 A and 605 1 A (see Royer et al. 1998 and their 
Table 1 for further details). The exposure times varied between 
1 and 2 minutes depending on the atmospheric conditions. 

The data were reduced independently by two of us (JM and 
YD) with different software tools. The results of both reductions 
were found to be in good agreement. The individual images were 
bias subtracted, divided by a median (calculated for each night 
and each filter), normalized flat-field and recentered in a ho- 
mogeneous way. We used the daophot software (Stetson '1987]l 
with aperture radii of 5, 7, 10, and 16 pixels (corresponding to 
3.5, 4.8, 6.9, and 11 arcsec). Whilst the results for the different 
apertures are generally in good agreement, we focus here on the 
light curves obtained for a 7-pixel radius aperture, which is well 
adapted to the general seeing conditions at the Haute-Provence 
site. The astrometric positions of the sources were obtained from 
combinations of the highest quality individual frames. Absolute 
(all-sky) photometry was derived from the larger aperture data, 
using a multi-night, multi-filter algorithm and a few standard 
stars (Manfroid I1993I I. This procedure yielded additional ref- 
erence stars for each field. These secondary standards and all 
non-variable stars were used to determine by a global minimiza- 
tion procedure, the zeropoints for the individual frames and each 
aperture radius, thus performing global differential photometry. 
Compared with the best-fit of the light curve (see below), we 
found a strange systematic trend in the O - C residuals of all 
filters, with the star progressively becoming fainter (by about 
0.02 mag) towards the end of the night (i.e. at higher airmasses). 
We carefully checked all reduction steps and found no reason 
for this behaviour No comparison star was found to exhibit the 
same behaviour In a Fourier periodogram of the residuals, we 
detected a peak at 1.01 d '. Whilst we were unable to identify 
the origin of this effect, which could thus be real, we note that 
(perhaps by coincidence) the effect disappeared when we dis- 
carded data taken at airmasses higher than 1.3. This is therefore 
the strategy that we adopt throughout this paper. 

2.2. Spectroscopy 

In September 1996 and April 1997, we obtained eight spectra 
of CygOB2 #5 with the echelle spectrograph mounted on the 
2.1m telescope at the Observatorio Astronomico Nacional of 
San Pedro Martir (SPM) in Mexico. The spectrograph covers the 
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spectral domain between about 3800 and 6800 A with a resolv- 
ing power of 18 000 at 5000 A. The detector was a SITe CCD 
with 1024 X 1024 pixels of 24finp-. The data were reduced us- 
ing the echelle context of the MroAS software, and specific orders 
covering important lines (Hj6, Ha, He i A 5876, He ii A 4686, C iii 
A 5696) were normalized using carefully chosen continuum win- 
dows. 

2.3. X-ray observations 

We obtained six XMM-Newton observations centred on 
Cyg OB2 #8 a, including Cyg O B2 #5 in the field of vie w of 
the EPIC cameras (Turner et al. SOOT Struder et al.'200P). The 
first four pointings separated by 10 days each were obtained in 
October - November 2004. Two follow-up observations were 
gathered in April - May 2007. The raw data were processed with 
the EPIC pipeline chains of the Science Analysis System (SAS) 
software version 6.0. Some bad time intervals characterised by 
high background events (so-called soft-proton flares) were re- 
jected. For a detailed description of the data reduction procedure, 
we refer to De Becker et al. (.2006i ). A few stray-light features 
(due to singly reflected photons) from CygX-3 are visible in the 
EPIC field of view. However, they do not affect the position of 
our target. We extracted the EPIC spectra of Cyg 0B2 #5 when- 
ever the source was unaffected by CCD gaps or bad columns. All 
the spectra were subsequently analysed using the xspec software 
(Arnaud fT996b . 

We also retrieved archival ROSAT-?S?C and A5CA-SIS ob- 
servations. Two PSPC data sets are available, rp200109n00 and 
rp9003 14n00. The latter was spread over a bit less than five days 
and we thus divided it into four time bins to search for orbital 
variability in the X-ray flux of Cyg OB2 #5. For A5CA, a single 
observation, obtained during the performance verification phase 
was obtained. This data set was already discussed by Kitamoto 
& Mukai ( I1996I I. All the archive data were reduced using the xs- 
ELECT package. Spectra of the source and the background were 
extracted, the latter over a source free region close to our target. 
The response matrices and ancillary response files were either 
retrieved from the archive or generated by the appropriate ftools 
routines. 

3. Results 

3.1. The optical iigint curve 

3.1 .1 . The ephemeris of Cyg 0B2 #5 

In the literature, the ephemeris of the primary eclipse was given 
by HJD = 2440 413.796 + 6.5977915 £ (Hall 1974) . Hubscher 
& Walter ( I2007I I reported the time of the minimum of the light 
curve as HJD 2453985.4928 + 0.0015 based on 15 measurements 
obtained in 2006. The time of minimum published by Hubscher 
& Wal ter ( |20071 l is 2057.006 cycles after the zero time of Hall 
( fT974] l. From our data, we infer HJD 2451049.702 + 0.025 to 
be the time of primary minimum, which corresponds to a shift 
of 0.04 in phase compared with Hall's ephemeris (1612 cy- 
cles). This apparent discrepancy prompted us to re-examine 
the question of the system's orbital period. Using the time of 
conjunction from the radial velocity curve of Wilson & Abt 
( |1951| l and the times of photometric minima listed by Miczaika 
( fT953T l. Hall ( fT974] l. Hubscher & Walter ( 120071 as well as 
our own data, we found that the best fit linear ephemeris was 
HJD = (2440413.750 ± 0.011) + (6.597840 + 0.000009) £. 
However, these new ephemerides lead to uncomfortably large 



phase shifts when applied to the epochs of our observations. 
Therefore, since the Hubscher & Walter (2007) result is ap- 
parently based on only 15 data points, we discarded this point 
from the fi0 A revised linear ephemeris would then become 
HJD = (2440413.718 + 0.013) + (6.597888 + 0.00001 1) E. At 
first sight these new ephemerides might look less precise than 
the old ones, but we emphasize that Hall ( 1974) did not provide 
an error bar in his ephemeris. Since the system is likely to be in 
an evolutionary phase with complex mass exchanges and mass 
loss (see below), we also considered the possibility of quadratic 
ephemerides. Using the times of photometric minima apart from 
the Hubscher & Walter (i2007i l value, we obtain 

1 . , 
HJD = HJDo +PqE + -P Po 

^ (2440413.796 + 0.024) + (6.597858 + 0.000014) E 
+(6.11 + 1.64) X 10"** £^ 

These are the ephemerides that we use throughout this paper 
The orbital period at the time of our photometric observations 
would hence be P = 6.598055 days, whilst the time derivative 
of the period would amount to P = (18.5 + 5.0) x 10"^ s s"' = 
(0.58 + 0.16) syr-i. 

3.1.2. The continuum filters 

We first studied the light curve of CygOB2 #5 in the ci and C2 
contmuum filterfl Restricting ourselves to observations taken 
at an airmass below 1.3, the total number of usable data points 
was 320 and 299 for the ci and C2 filters, respectively. The data 
provided a good coverage of the orbital cycle, although they did 
not sample the core phase of secondary minimum. These light 
curves indeed confirm that both minima have different depths 
and also indicate a different brightness immediately after the 
primary eclipse compared to the situation after the secondary 
eclipse. These features agree with the description of Hall ( 1974|l 
and cannot be explained by stars with a uniform surface bright- 
ness. We thus considered situations where at least one of the two 
stars has a non-uniform brightness distribution (see below). To 
analyse this light curve further, we used the nightfall softwar^ 
developed and maintained by R. Wichmann, M. Kuster and P. 
Risse. 

Rauw et al. (1 19991 1 derived two slightly different values of the 
dynamic mass ratio q - mp/nis - 3.31 ±0.25 when only the new 
RV data were used in the orbital solution and q = 3.54 + 0.18 
when literature RV data were also included. In a first series of 
trials, we hence tested the capabilities of the light curve analysis 
to distinguish between these different values of the mass ratio. 
We found that within a 3 cr interval around the mean mass ra- 
tios, the light curve fits were not very sensitive to the value of q. 
However, the fits for the lower value of ^ = 3.31 are of some- 
what better quality and in the following we thus keep the mass 
ratio fixed at the value derived from the new RV data of Rauw et 
al. (1999). 

At first, we tested models with the surface temperature dis- 
tribution ruled by gravity darkening only (i.e. no bright or dark 
spots). The NIGHTFALL code immediately converges towards a 

- Given the difficulties in defining accurately the time of primary 
minimum in our data (containing several hundred measurements), the 
error in the time of minimum quoted by Hubscher & Walter ( |2007t 
seems extremely optimistic to us. 

^ Table 1 with the light curves of Cyg 0B2 #5 in the various filters is 
available in electronic form at the CDS. 

* ,http://www.hs.uni-hamburg.de/DE/Ins/per/Wichmann/Nightfall.html| 
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contact (or even over-contact) configuration where both stars fill 
up their Roche lobe and are seen under an inclination ; close to 
62 - 67°. The best-fit model temperature of the secondary star is 
29 500 K in this case (assuming a primary effective temperature 
of 36 000 K). 

The residuals of this solution indicate that the data points 
during the egress of the primary eclipse are systematically lower 
than the model predictions, whilst the agreement between the 
data and the model is reasonable over the remainder of the light 
curve. We note that allowing for a slightly eccentric orbit does 
not solve this discrepancy between the data and the synthetic 
light curve. This feature could be explained instead by either the 
secondary being somewhat brighter on its leading side or the 
primary being fainter on its leading side. Since the first option 
is reminiscent of the location of the colliding wind region dis- 
cussed by Rauw et al. ( 1999 ), we tested this model by adding a 
bright "spot" on the secondarjQ This bright spot model yields a 
significantly better fit. However, there is a degeneracy between 
the secondary surface temperature, spot size, and spot tempera- 
ture ratio. These parameters combine in such a way as to produce 
a mean temperature of the secondary star over its entire surface 
of about 35 600 K. As a result, the secondary's surface temper- 
ature outside the spot area is no longer constrained and we thus 
decided to adopt the value of 29 500 K derived above. We then 
find that the part of the secondary star facing the primary is sig- 
nificantly hotter (the hot spot has a best-fit model temperature of 
36 500 K, as compared to 29 500 K for the rest of the secondary's 
surface) and hence brighter (by a factor of 2.35 in bolometric lu- 
minosity). Actually, the surface of the "spot" covers more than 
half of the total surface of the secondary star! However, the cen- 
tre of this hot spot is found not to be aligned with the binary axis. 
In fact, the best-fit model yields a longitude for the spot centre 
of -15°, where 0° corresponds to the binary axis and negative 
longitudes refer to the leading side of the secondary's surface. 
The fact that the centre of the spot is located on the leading side 
of the secondary star suggests indeed that it cannot be due to 
a sole heating effect by the primary's radiation, but that addi- 
tional heating from a Coriolis-deflected wind interaction region 
wrapped around the secondary star is probably playing a major 
role. This situation is in excellent qualitative agreement with the 
expectations for the wind interaction zone characterised by vari- 
ations in the equivalent width of the He ii A 4686 emission (Rauw 
et al. lT999l l. 

The mean values of the absolute differences between the ob- 
served and computed magnitudes \0 - C\ are 0.009 mag for both 
filters. These values are slightly, but significantly larger than the 
standard deviations (rms) of the differential magnitudes of non- 
variable stars of comparable magnitude, probably revealing in- 
trinsic variability of the system in addition to the systematic or- 
bital effects. 

To establish the error bars on the best-fit parameters in 
Table|2] we explored the parameter space by keeping one of the 
parameters fixed and allowing the others to vary until a mini- 
mum in was reached. The error bars were then calculated for 
the variation in corresponding to a 90% confidence level and 
the appropriate number of degrees of freedom. We emphasize 
that the parameter space is complex and some parameters (such 
as the spot's temperature ratio and its size) are correlated to some 
extent. 



^ We note that we assumed the spot to be on the stellar equator, which 
should be a reasonable assumption since the hot region is likely to arise 
from mutual heating effects and/or interactions with the primary's wind. 



Table 2. Best-fit parameters of the continuum light curves of 
CygOB2#5. 



Parameter ci & C2 

7n 64.2!l-^ 

q 3.31 (fixed) 

Filling factor primary 1.00^ Jj| 

Filling factor secondary 1.00^ Jjj 



r,ffp(K) 36000 (fixed) 

reff.s (K) 29500 (fixedj 



Bright "spot" 




Longitude 


(-15!!^)= 


Latitude 


0° (fixed) 


Radius 


(130!i^)° 


Temperature ratio 


1.24f?i 



" This value was derived from the best-fit model without a spot on 
the secondary star and refers here only to the un-spotted part of the 
secondary's surface. 




(/) (P - 6.598055 days) 

Fig. 1. The ci and C2 light curves of Cyg OB2 #5 along with 
our best-fit model (Table|2]l. The phases were computed using 
our quadratic ephemeris, and phase = 0.0 corresponds to the 
primary minimum (secondary star in front). The panels below 
the light curves illustrate the deviations O - C for each filter. 



From the orbital solution of Rauw et al. (1 19991 1 and the in- 
clination of 64.2°, we infer the absolute masses of the compo- 
nents of CygOB2 #5 to be (31.9 + 3.2) Mq for the primary and 
(9.6 + 1.1) Mq for the secondary. With mean (averaged over the 
entire surface) surface temperatures of 36 000 and 35 600 K for 
the primary and secondary respectively, the bolometric luminos- 
ity ratio between the two stars is about 3.1, whilst the optical 
(V-band) brightness ratio should be around 2.9. The total bolo- 
metric luminosity of the system equals about log Lboi/io = 6.04 
corresponding to an absolute magnitude of My - -7.05. The 
observed V magnitude outside eclipse is 9.20 around maximum 
(Hall [T974] l, whilst it is 9.80 for the ci filtefl For the reddening 
of CygOB2 #5, Torres-Dodgen et al. (fTWD l quoted Ay = 6.40, 



* Using the Cardelli et al. ( |1989t reddening law, we infer A^i/Ay = 
1.106. 
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Fig. 2. The light curves of Cyg OB2 #5 as observed through the 
Une-bearing filters. The solid line yields the best fit solution of 
the ci and C2 filters (Table|2|. The panels below the light curves 
illustrate the deviations between the actual data and the best fit 
solution for the continuum filters. 



whilst Wegner (120031 ) gave a value of Ay - 6.35. These numbers 
also agree with the average B -V - 1 .68 value of the binary in- 
ferred from the data compiled by Reed (I1998I I. Hence, our anal- 
ysis of the light curve yields a distance modulus in the range 
9.77 - 9.90 corresponding to a distance of about 900 - 950 pc. 
This is significantly less than the value of Torres-Dodgen et al. 
( 1199 1'). who used ubvyJHKL' photometry of the Cyg OB2 asso- 
ciation to infer a distance modulus of 11 .2 + 0.2. The latter value 
is in good agreement with the distance modulus of 1 1 .2 proposed 
by Massey & Thompon (1991 ). It is obvious that our result de- 
pends strongly on the (assumed) temperature of the primary and 
secondary components. Here, we adopted 36 000 K for the pri- 
mary as inferred by Rauw et al. ( 1999 ) from the 06.5-7 spectral 
type. However, to ensure that our results agree with the 1.7kpc 
distance estimate, we would require a much hotter primary (of 
order 48 000 K) which appears extremely unlikely. On the other 
hand, our distance estimate is not very sensitive to uncertainties 
in the secondary's temperature. 

Whilst our best-fit parameters differ significantly from those 
of Leung & Schneider ( 1978 ), we emphasize that a similar dis- 
crepancy between the luminosity of the light curve fit and that 
inferred from the typical distance was also noted by the latter 
authors, although they formulated this problem in a different 
way. We therefore encourage other studies of eclipsing binaries 
in Cyg OB2 (e.g. Cyg OB2 #3; see Kiminki et al. l2008l l to con- 
firm independently the distance of the cluster. 

3.1.3. Line-bearing filters 

We then compared the best-fit solutions of the continuum filters 
with the observed light curves of the filters centred on the lines. 
The results are shown in Fig.|2] For each of the three line-bearing 
filters, we simply shifted the theoretical c\ or C2 light curve un- 



til it fitted the average magnitude measured in the line-bearing 
filter The remaining residuals then inform us about effects that 
are specific to a given filter No strong effect is evident for filters 
centred on Hen /1 4686 (306 data points, \0 - C\ = 0.008 mag) 
and Civ ^^5801, 5812 (44 data points, \0-C\ = O.OlOmag). 
The average residuals are fully consistent with the results ob- 
tained for the c\ and filters. In fact, the spectrum of Cyg OB2 
#5 does not contain any strong emission or absorption lines in 
the Civ /I/15801, 5812 filter. Hence, it is unsurprising that this 
filter behaves in the same way as the continuum filters. However, 
a very different situation is observed for the filter centred on 
He I ^5876, for which we obtain \0-C\ = 0.022 mag (312 data 
points) and Fig.|2]clearly reveals a phase dependent modulation 
of the residuals. Actually, the light deficit after primary mini- 
mum appears stronger in this filter than in any other filter that 
we have used. As a result, the residuals indicate a broad mini- 
mum around phase 0.15. We therefore attempted to fit the light 
curves of the three line-bearing filters with the same model as 
for the continuum filters, but allowing the secondary's tempera- 
ture as well as the spot's temperature ratio, longitude, and size to 
vary until we achieve the best-fit for each filter The results are 
listed in Table[3] 

Table 3. Parameters of the best-fit solutions to the light curves 
of Cyg OB2#5 as observed through the line-bearing filters. 



Filter" 


Hen 


Civ 


He I 


Teff,., (K) 


28400 


32600 


21300 


Bright spot 








Longitude (°) 


-15 


-28 


-26 


Radius (°) 


125 


117 


118 


Temperature ratio 


1.14 


1.26 


1.46 



" Several parameters were frozen in these fits: ; 
primary and secondary filling factors = 1.00, Teff p 
of the bright spot = 0°. 



= 64.2°, q 
36000 K, 



= 3.31, 
latitude 



From this table, we find that the apparent bolometric bright- 
ness ratio of the spot is essentially the same in the continuum 
(2.36), He II (1.69) and Civ (2.52) filters, but reaches 4.54 in the 
He I filter. We also note that the data for the latter filter yield the 
lowest secondary temperature, i.e. in this filter the primary ap- 
pears much brighter than the unspotted part of the secondary's 
surface. We propose an explanation for this behaviour in the next 
section. 



3.2. Optical spectroscopy 

Unfortunately, our optical spectra do not provide a complete 
coverage of the orbital cycle. Nevertheless, we can use them to 
illustrate the variability of the most prominent lines (see Fig.|3]l 
and attempt to understand the origin of the discrepant behaviour 
of the He i A 5876 filter in our photometric data. 

The variations in the He ii A 4686 and Ha emissions are very 
similar to the variability pattern reported by Rauw et al. d 1 9991) 
and Vreux ( 1985 ), respectively. We mapped the emission regions 
of these emission lines in velocity space using the Doppler to- 
mography technique of Rauw et al. (2002) . For this purpose, 
we adopted a reference frame centred on the centre of mass of 
the binary with the x-axis pointing from the primary to the sec- 
ondary and the positive y-axis pointing along the direction of the 
secondary's orbital motion. The basic assumption of a Doppler 
tomography technique is that the phase dependence of the radial 
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Fig. 3. Line profile variability in the optical spectrum of Cyg OB2 #5. The panels illustrate from left to right the observed variations 
in the He ii /1 4686, He i A 5876, and Ha lines. The successive spectra are shifted vertically by 0.2 (respectively 0.3) continuum units 
in the left and middle (respectively the right) panels. In all cases, the profiles are shown in the heliocentric frame of reference and the 
intensities of the lines were corrected for the variations in the continuum light curve. In the left and middle panels, the short-dashed 
lines show the central wavelength of the narrow-band filter centred on these lines, whilst the two long-dashed lines provide the 
FWHM of these filters. The phases of the observations are indicated to the right of the figures (quadratic ephemeris with = 0.0 
corresponding to the primary minimum). 



velocity v((p) of any gas flow that is stationary in the rotating 
frame of reference of the binary can be described by a so-called 
'S-wave' relation: 

v(<p) - -Vx COS {2n(p) + Vy sin {2n(p) + v^ (1) 

where (p is the orbital phase (<^ = at primary eclipse), whilst 
Vx and Vy are the projected velocity components along the x and 
y axes and Vj corresponds to the apparent systemic velocity of 
the line under investigation (assumed to equal the primary's sys- 
temic velocity of -33kms"'). The data were weighted to ac- 
count for the uneven sampling of the orbital cycle, and the inten- 
sities of the emission lines were corrected for variations in the 
underlying continuum. 

While the limited sampling of our data leads to some low- 
level artefacts in the Doppler map (suppressed in Fig. |4] by 
choosing the lowest contour level at 0.50 x the intensity of the 
highest peak), it can be seen that there is good qualitative agree- 
ment between the results obtained here and the simple S-wave 
analysis of Rauw et al. (119991 1. For the Hen /i4686 line, we 
found a single emission region concentrated close to the veloc- 
ity of the centre of mass of the secondary star, although with 
a slightly positive v^. The Doppler map of the Ha line is far 
more complex, reflecting the multi-peak structure of the profile 
of this line. Nevertheless, the bulk of the emission originates in 
material with positive and Vy velocities. These emission lines 
are thus formed mostly in material that moves from the primary 
towards the secondary and in the same orbital direction as the 
secondary. Although more complex configurations cannot be ex- 
cluded at this stage, we note that this picture agrees with a col- 
liding wind scenario where the primary's wind dominates the 
secondary wind. In such a case, the optical emission lines would 
arise from recombination in higher density regions due to the 
impact of the primary wind on the secondary's wind or photo- 
sphere. This higher density region would be geometrically coin- 
cident with the brighter and hotter region of the secondary star 
found in our photometric analysis. This description is similar to 
that presented by Rauw et al. (1999), although we caution that 
these authors assumed a detached configuration, which is clearly 
not the case here. The situation in Cyg OB2 #5 is thus quite dif- 
ferent from a 'classical' wind - wind collision. Indeed, according 



to our photometric solution, the stars are in contact and the wind 
interaction mainly concerns those parts of the secondary surface 
that sweep up the probably more energetic primary wind during 
the orbital motion. Given the proximity between the components 
of the binary system, the winds very probably interact at low ve- 
locities and the shock-heated material is thus unlikely to reach 
temperatures sufficiently high to produce a sizeable X-ray emis- 
sion. 

The He i A 5876 line displays complex variations going from 
a P-Cygni type profile at phases between 0. 1 and 0.4 to an almost 
pure emission line over the second half of the orbital cycle (see 
Fig-Ell. This explains the strange behaviour of the photometric 
data from the filter centred on this line as discussed in the previ- 
ous section. In fact, the equivalent width of the line, corrected for 
the photometric variations in the continuum, changes from about 
1.0 A (for the absorption-dominated P-Cygni type profiles) to 
-1.2 A for the emission-dominated profiles. This is quite a sig- 
nificant variation, sufficient to produce the 0.08 mag (peak-to- 
peak) deviation of the Hei /1 5876 filter from the pure contin- 
uum filters. The emission component of the line hence probably 
arises from a region of limited angular extension, most likely 
over or above the hot spot on the secondary star Unfortunately, 
the highly variable absorption component in this profile prevents 
us from applying the Doppler tomography technique to the He i 
^5876 line. 

In Fig. [3] we note some hints of slight differences between 
the profiles observed in 1996 and 1997 at similar orbital phases. 
In 1997, the absorption trough of the P-Cygni profiles displays a 
rather sharp absorption component at phase ^ - 0.13, whilst no 
such feature is seen at - 0.10 in the 1996 data. This could po- 
tentially be related to long-term variations similar to those seen 
in the Hp line profile variations (Rauw et al. |1999l l. 

3.3. X-ray data 

We analysed all the EPIC spectra of CygOB2 #5 using the 
xsPEC software. The spectra were fitted with a two-temperature, 
optically-thin, thermal plasma model mekal (Mewe et al. 119851 
Kaastra [1992) absorbed by a fixed, neutral interstellar absorp- 
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Fig. 4. Doppler maps of the He n /1 4686 (top panel) and Ha (bot- 
tom panel) lines in the spectrum of Cyg OB2 #5. The crosses cor- 
respond to the radial velocity amplitudes of the centre of mass 
of the primary and secondary. The shape of the Roche lobe in 
velocity-space (thick dashed line) has been calculated for a mass 
ratio (primary/secondary) of 3.31. The black dot indicates the 
position of the highest peak and the contours correspond to lev- 
els of 0.95, 0.80, 0.65, and 0.50 times the maximum emissivity. 



tion column as well as a column of (partially) ionized circum- 
stellar wind material. Within xspec, we hence used a model of 
the type wabsxwindxCmekal + mekal). The best-fit parame- 
ters for the XMM-Newton-EPlC spectra are listed in Tableland 
illustrated in Fig.|6]for models that assume a solar composition. 
The first column of Table|4] yields the dates of the observations 
(at mid exposure) with the uncertainty corresponding to half the 
integration time. The interstellar column-density of neutral hy- 
drogen was fixed to the value 1.13 x 10^^ cm"^ derived from the 
E(B - V) - 1 .94 colour excess inferred by Torres-Dodgen et al. 
( I1991I I and the mean relation between A^h and E(B-V) of Bohlin 



et al. ( I1978I I. The wind component of the absorption was mod- 
elled using the opacities computed with the model developed by 
Naze et al. (2004). The norm of the mekal model is defined to be 
'°4;rrf™ where d is the distance (cm) and EM - J rieriH dV is 
the emission measure (in cm"^). The last two columns of Table|4] 
yield the observed and ISM absorption corrected fluxes in the 0.4 
- 10.0 keV energy domain. 



> 



0,1 ^ 



0.01 E 



0.001 




Energy (keV) 

Fig. 5. The EPIC-pn (black), -MOSl (red) and -M0S2 (blue) 
spectra of CygOB2 #5 as observed on JD 2453328.543 (18 
November 2004). The data points are shown with their error 
bars, whilst the continuous line yields the best-fit model (see 
Table|4|i. The bottom panel yields the contributions of the in- 
dividual energy bins to the of the fit with the sign of the cor- 
responding O - C difference. 



We tested single temperature fits, but these failed to repro- 
duce the hard part of the spectrum. Significantly better results 
were achieved for two-temperature fits (see Fig.|5]l. The lower 
temperature is found consistently to be around 0.63 ± 0.01 keV, 
whilst the hotter component at kT2 = (1.74 ± 0.14)keV is some- 
what less constrained. The observed flux varies by more than a 
factor of two. This is unlikely to be due to a modulation of the 
circumstellar absorption, since (1) N^ind varies relatively little 
and (2) N„ind is at its second highest value when the observed 
flux reaches its maximum. According to our fits, the flux varia- 
tions are most likely due to a changing emission measure of the 
X-ray emitting plasma. 

We compared the X-ray fluxes of the various pointings over 
seven narrow energy bands (0.25 - 0.6, 0.6 - 1.0, 1.0 - 1.5, 1.5 
- 2.0, 2.0 - 3.0, 3.0 - 5.0, and 5.0 - lO.OkeV). Whilst the four 
XMM-Newton spectra taken in 2004 do not show variations of 
more than 10% in any of the energy bands, the fluxes appear 
much higher in the two 2007 spectra, in some bands by more 
than a factor of 2. Comparing the last two observations with the 
first four, the increase is by far the most important at energies 
above 1 keV and increases with energy. This underlines that the 
spectrum in our 2007 pointings is significantly harder than those 
in 2004 as can already be seen by the increase in the norm of the 
hotter component in the two-temperature fits in Tabled 

To answer the question of whether or not the observed vari- 
ability is related to the orbital period of the binary system, we 
also analysed archival ROSAT and ASCA spectra. 
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Table 4. Results of the best-fit models of the EPIC spectra of CygOB2 #5. 



JD - 2400000 logN„ind kT\ norm; kT^ norm^ xt (d-o.f.) f^' ff^ 

keV keV ergcm^^s"' erg cm"- s 



53308.579 


±0.122 


21 53+°'^ 


0.64! 


02 
02 


(9.04!i;g) X 10-3 


1.66!H 


(2.62! 32) X 10- 


-3 


1.40 (472) 


2.66 X 10- 


12 


1.46 X 10- 


ll 


53318.558 


±0.133 


21.63!;;^ 


0.61! 


02 
02 


(12.04!|^^)x 10-3 


1.64!;^ 


(2.93! ^^) X 10- 


-3 


1.24 (432) 


2.98 X 10- 


12 


1.60 X 10- 


11 


53328.543 ±0.145 


21.56+-^^ 


0.63! 


02 
01 


(10.89!»|^)xlO-3 


1.98+|» 


(2.19!^6)x 10- 


-3 


1.06 (580) 


2.84 X 10- 


12 


1.58 X 10- 


11 


53338.506 


±0.133 




0.62! 


02 
02 


(9.68!];^^) X 10-3 


i.59!;i^ 


(3.23!^^^) X 10- 


-3 


1.19 (384) 


2.91 X 10- 


12 


1.60 X 10- 


11 


54220.366 


±0.167 


21.79!«J 


0.65! 


03 
02 


(16.23!^^|)x 10-3 


1.76!?^ 


(4.99!^^) X 10- 


-3 


0.91 (458) 


4.37 X 10- 


12 


1.69 X 10- 


11 


54224.180 


±0.173 


21.76!;!!* 


0.63! 


02 
02 


(19.99!|^^) X 10-3 


1.79!]2 


(6.39!-^g) X 10- 


-3 


1.06 (614) 


5.51 X 10- 


12 


2.20 X lo- 
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Fig. 6. The main parameters of the fits to the EPIC data are 
shown as a function of phase (quadratic ephemeris with (p - 0.0 
corresponding to the primary minimum). From top to bottom, 
the panels illustrate the observed flux in the 0.4 - lO.OkeV band 
(in units lO-'^ergcm-^s-'), the circumstellar (wind) column 
density and the two temperatures (kT\ and kT2) respectively. 



Previous studies of these data yielded conflicting conclu- 
sions about the X-ray variability of CygOB2 #5. Comparison 
between the ROSAT observations and EINSTEIN data in fact 
suggested that the star had apparently shown two events of en- 
hanced (increase by a factor 2 in flux) X-ray emission between 
late 1979 and early 1980 (see Waldron et al. 119981) as well as 
in the mid 1990s (ROSAT All Sky Survey data). Waldron et al. 
( I1998I I pointed out that since the star was in its low level ra- 
dio state during the first event and in a high level state during 
the second event, these X-ray events are unlikely to be related to 
the radio emission. Finally, Corcoran (1996) analysed the phase- 
dependence of the soft X-ray count rate from pointed ROSAT- 
PSPC observations. No significant modulation was found for 
Cyg OB2#5. 

We thus analysed flie ROSAT-PSPC and A5CA-SIS spectra, 
using a single temperature model for the former and a two- 
temperature model for the lattefl These fits were then used to 



' Given the limited energy range of the PSPC instrument, a second 
temperature would not be constrained by the data and is actually not 
required to fit the spectrum. However, the hard part of the SIS spectrum 



evaluate the observed fluxes in the 0.5 - 2.0 keV domain. The re- 
sults are listed in Table|5] The first and second columns yield the 
dates of the observations (at mid-exposure and with the uncer- 
tainties corresponding to half the integration time) and the corre- 
sponding orbital phases computed from our quadratic ephemeris. 
The third column indicates the instrument used in the observa- 
tion and the last column provides the observed flux in the 0.5 - 
2.0 keV energy range. 

As far as the ROSAT-PSPC data are concerned, we find 
indeed surprisingly little variability: the flux remains constant 
within 5% of its mean level. At orbital phases that were observed 
with all three instruments, the agreement between PSPC, SIS, 
and EPIC fluxes was extremely good given the very different 
sensitivities of these devices and the fact that the A5CA-SIS ob- 
servation is limited by stray-light contamination from Cyg X-fl 



Table 5. The soft (0.5 - 
function of time. 



2.0 keV) X-ray flux of Cyg OB2 #5 as a 



JD - 2400000 







Instrument 


robs 

erg cm-2 s 


-1 


48368.064 


±0.137 


0.570 


±0.021 


PSPC 


1.59 X 10- 


12 


49106.963 


± 0.466 


0.558 


±0.071 


SISO& SISl 


1.66 X 10- 


12 


49107.121 


± 0.087 


0.582 


±0.013 


PSPC 


1.48 X 10- 


12 


49109.234 


± 0.058 


0.902 


± 0.009 


PSPC 


1.63 X 10- 


12 


49110.131 


± 0.203 


0.038 


±0.031 


PSPC 


1.62X 10- 


12 


49111.230 


± 0.376 


0.205 


± 0.057 


PSPC 


1.64 X 10- 


12 


53308.579 


±0.122 


0.354 


±0.018 


EPIC 


1.39 X 10- 


12 


53318.558 


±0.133 


0.866 


± 0.020 


EPIC 


1.57 X 10- 


12 


53328.543 


±0.145 


0.380 


± 0.022 


EPIC 


1.46 X 10- 


12 


53338.506 


±0.133 


0.890 


± 0.020 


EPIC 


1.51 X 10- 


12 


54220.366 


±0.167 


0.543 


± 0.025 


EPIC 


1.98 X 10- 


12 


54224.180 


±0.173 


0.121 


± 0.026 


EPIC 


2.44 X 10- 


12 



From Fig.|2l we note that the XMM-Newton observation with 
the highest X-ray flux was taken when the secondary star was 
moving away from us. However, it seems rather unlikely that we 
are dealing with a phase-locked variability since the EINSTEIN- 
IPC observation in which the system was detected at a similarly 
high flux level (Waldron et al. 1 19981 1 was spread over five days 
whilst the ROSAT All Sky Survey data (when the star was also 



can only be fitted by including a second temperature, although the latter 
is poorly constrained due to the limited quality of the SIS data. 

We emphasize that, to within 5%, we obtain the same flux as found 
by Kitamoto & Mukai ( fT996t in their analysis of the A5CA-SIS CCD 
spectra. 
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Fig. 7. The soft band (0.5 - 2.0 keV) flux of CygOB2 #5 as a 
function of orbital phase of the eclipsing binary (top panel) and 
as a function of time (bottom panel). The filled circles indicate 
OUT XMM-Newton EPIC data, the triangles stand for the ROSAT- 
PSPC data whilst the open square corresponds to the AS'CA-SIS 
observation. The dashed lines indicate the epochs of the mini- 
mum radio flux at 6 cm according to Kennedy et al. (120091 ) . 



at a high X-ray level, Waldron et al. '1998") were collected over 
an even longer time interval. Therefore, we conclude that there 
is currently no evidence for a phase-locked modulation of the 
X-ray flux, at least not with the period of the eclipsing binary 
system. 

An interesting result is the lack of eclipses in the X-ray do- 
main. Indeed, our observations reveal no change in the observed 
X-ray flux near the phases of the eclipses (Fig.|2ll. Therefore, the 
X-ray emitting region must be far wider than the typical dimen- 
sions of the stars in the eclipsing binary. A similar situation was 
reported for jSLyrae (B6-8IIp -i- BOV, Ignace et al. 2008), for 
which it was suggested that the X-rays arise in the wind of the BO 
mass-gainer star and that the constancy of the X-ray flux is due 
to scattering in an extended halo from above and below the disk 
around the BO star (Ignace et al. 2008 ). In the case of Cyg OB2 
#5, the X-ray emission could originate in various places. First, 
it could be intrinsic to the outer regions (several stellar radii) of 
the primary star's atmosphere. However, the clear detection of a 
rather hard component in the spectrum argues against this expla- 
nation. It could also, at least partially, originate from the wind- 
wind interaction between the two stars. Since the stars are in a 
contact configuration, we emphasize that we would expect this 
X-ray emission to be rather soft if anything, because the winds 
would interact at velocities far lower than the terminal velocity. 
We would therefore expect at least some modulation of the wind 
column density between phases 0.0 (secondary wind in front) 
and 0.5 (primary wind in front). The existence of a hard compo- 
nent and the lack of significant modulation of the wind column 
density thus argue against this scenario. The most likely inter- 
pretation is then that at least the hard part of the X-ray emission 
originates in the interaction region between the combined stellar 
winds of the eclipsing binary and the wind of the tertiary compo- 



nent or the wind of a yet unknown fourth component, responsible 
for the 6.7 yr modulation of the radio emission (Kennedy et al. 
120091 ). As far as the former possibility is concerned, it remains 
to be verified whether the emission measure of this interaction 
region in such a wide system with a period of centuries would 
be sufficient to produce a significant hard component in the spec- 
trum. Concerning the latter possibility, the plot of the soft X-ray 
flux as a function of time in Fig.|7] shows that during the 6 cm 
radio minimum state (depicted by the dashed lines in Fig.|7]), we 
observe both, low and high X-ray fluxes. Therefore, it seems un- 
likely a priori that the X-ray variations have the same origin and 
timescale as the radio variations. We also note the striking simi- 
larities between the X-ray spectra of Cyg OB2 #5 and the triple 
system HD 167971 that consists of a close eclipsing binary with 
a third component with a period of several years (see De Becker 
et al. 12005] ). 

If we adopt a distance of 925 pc, the average (ISM-corrected) 
X-ray luminosity in the 0.4 - 10.0 keV band would be 1.73 x 
lO-'^ergs Compared to our estimate of the bolometric lumi- 



-6.39. 



nosity of the eclipsing binary, we thus obtain log 
The system thus appears slightly (by a factor of 3) X-ray overlu- 
minous compared with the Lx-L\,o\ relation established by Sana 
et al. ('2006 ) for the O-type stars in NGC 623 1 . Whilst this is only 
a mild overluminosity, it is consistent with our interpretation that 
part of the X-rays originate from wind - wind interactions. 

Finally, we have also used the XMM-Newton data to search 
for short-term X-ray variability (on the timescale of individ- 
ual pointings). Hydrodynamical instabilities in the plasma in 
a colliding-wind region could in principle create these short 
timescale variations (Stevens et al. 1992). We considered the 
EPIC light curves with binnings of 100, 500, 1000, and 4000 s, 
in three diff'erent energy bands: [250: 1000] eV, [1000:2500] eV, 
and [2500: 10000] eV. We applied three independent statistical 
tests to the data (x^, Kolmogorov-Smirnov and probability of 
variability test; see Sana et al. 120041 ) to detect variations on the 
timescale of one observation. No significant variability could be 
detected in any of the energy bands, for any instrument, and/or 
any time binning. 



4. Discussion 

4.1. Evidence for a changing period? 

We begin by considering the meaning of the non-linear 
ephemeris derived from the time of minima. First of all, we 
emphasize that a linear ephemeris definitely does not pro- 
vide a good fit to the times of primary minimum. However, 
since there are only a limited number of extensive (i.e. cov- 
ering several cycles) observing campaigns of this system, the 
quadratic ephemeris clearly requires confirmation. This is es- 
pecially true since V 382 Cyg, another eclipsing O + O over- 
contact binary, apparently displays some stochastic variations 
in its period (Harries et al. 1 19971 and references therein) and 
we cannot exclude the existence of a similar phenomenon in 
the case of our target. This being said, we have found that the 
time derivative of the orbital period of Cyg OB2 #5 amounts to 
P = (18.5 + 5.0) X lO^'^ss-' = (0.58 + 0.16)syr-'. Although 
this value appears rather large, it is comparable to those values 
inferred for other evolved eclipsing binaries such as V 444 Cyg 
(WN5 H- 06, P = 0.12 syr-^ Eaton & Henrv [19941 ) , GPC ep 
(WN6/WCE + 03-6, P = 1.3 syr '; Demers et al. "2007) or 
CQCep (WN7 + 07, P = -0.014 syr '; Antokhinaet al. 1987]). 

Wellstein (2001) pointed out that the detection of changes 
in the orbital period of CygOB2 #5 could help constrain the 
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mass-loss rate and hence the evolutionary stage of this system. 
However, the interpretation of the rate of period change de- 
pends strongly on the way that mass is exchanged inside or lost 
from the system. Depending on the details of the process, the 
orbital period can either increase or decrease as a function of 
time (KhaliulUn [19741 Singh & Chaubev fT986] l. If we assume 
that mass is transferred from the primary towards the secondary 
without any change in the orbital angular momentum and the 
systemic mass (case 1 of Singh & Chaubev I 19861 1. we find that 

1 mpMs P 
3 (mp - nis) P 

This then yields a mass-transfer rate of ihp - + 1.5) x 

10"^ Mo yr"', i.e. mass would be transferred from the secondary 
towards the primary star. As we argue below, this seems rather 
unlikely. 

Alternatively, we can assume that the mass is lost primarily 
through an isotropic (to first order) stellar wind (case IV of Singh 
& Chaubey 11986b . In this situation, the mass-loss rate can be 
directly computed to be 

1 P 
nip^-- {nip + m,) - 

In our case, we thus obtain a mass-loss rate from the binary 
system of nip - -(2.1 + 0.6) x lO^^Mgyr"' quite similar to 
the value inferred for the Wolf-Rayet component of GPCep 
(M - -2.9 10"^ Mq yr"' ; Demers et al. 120021 ) and about a factor 
of 3 higher than the most extreme values found for V 444 Cyg 
(Eaton & Henry lT994l l. 

We emphasize that in the particular case of CygOB2 #5 
an additional effect could alter the observed times of minimum 
light. This effect is related to the third component dicussed by 
Contreras et al. (1 19971 ). As described in Sect.[T] these authors 
confirmed the presence of a third star, probably of spectral type 
BO-2 V, at an angular distance of 0.98 arcsec from the close bi- 
nary. If this third star is gravitationally bound to the binary sys- 
tem and if the angular separation corresponds to the radius of 
a circular orbit or the semi-major axis of an elliptical orbit, the 
corresponding orbital period would be about 9200 yr for a dis- 
tance of 1.7 kpc (or 3700 yr for a distance of 925 pc). The sep- 
aration between the epoch of Hall (119741) and our data is about 
29.1 yr. If we make yet another assumption that the triple sys- 
tem is seen under an inclination of 90° and is currently close to 
quadrature phase, then the maximum change in the distance of 
the eclipsing binary towards the Earth would be about 5x10"^ km 
(or about 7x10^ km if the system lies at 925 pc). Due to the fi- 
nite velocity of light, such a change in distance can produce a 
lag in the time of minimum light. Interpreting the 0.266 day de- 
lay between Hall's linear ephemeris and our determination of 
the time of minimum in this way, requires a change in distance 
of 6.9 x 10"^ km. Whilst there are admittingly a large number of 
assumptions here, we nevertheless note that this light-travel ef- 
fect could contribute significantly to the observed deviation from 
the old linear ephemeris. 

In summary, it is probable that several effects (mass-transfer, 
mass-loss, and orbital motion in a triple system) conspire to pro- 
duce a non-linear ephemeris in Cyg OB2 #5. Additional accurate 
determinations of the primary minimum over a long timescale 
would be most welcome in clarifying this issue. 

4.2. The interactions in fhe eclipsing binary system 

An important result of our study is that the components of the 
eclipsing binary are currently in a contact configuration. In ad- 



dition, we find that the secondary appears hotter over part of its 
surface and this is probably due to the energy transfer in a com- 
mon envelope and the heating of the secondary star by both the 
primary's radiation field and the impact of the primary's stellar 
wind on the secondary's surface. It is interesting to compare this 
configuration with the radiation-hydrodynamics simulations of 
Dessart et al. (2003 ). These authors simulated mass transfer in a 
binary system in the presence of strong stellar winds and showed 
that neither the stellar wind's momentum nor the radiation field 
of the accreting star can prevent mass transfer by Roche lobe 
overflow. Dessart et al. (120031) further showed that the wind of 
the accreting star is always dominated by the primary's wind and 
the mass stream. The fact that the Doppler maps of Cyg OB2 #5 
and the S-wave analyses performed by Rauw et al. ( 119991 ) in- 
dicated that (line-emitting) material is moving towards the sec- 
ondary hence suggests that the primary is transferring matter to 
its companion. 

Our finding that the X-ray flux of the system is not modu- 
lated by the short period of the eclipsing binary indicates that the 
bulk of the X-ray emissions does not originate from the wind- 
wind interaction between the two components. It rather seems 
that the interaction between the wind of the eclipsing system and 
that of a tertiary or maybe even a fourth component could be at 
the origin of the X-ray emission. It is likely that this also holds 
for the acceleration of relativistic particles responsible for the 
non-thermal radio emission (Kennedy et al. 2009 ) and the poten- 
tial contribution to the yet unidentified EGRET y-my (lOOMeV 
-20GeV) source 3EG J2033+4118 (Benaglia et al. .200lfl . 

4.3. Tine evolutionary state 

An important consequence of our finding that the system is 
currently in a contact configuration (as originally suggested by 
Leung & Schneider 1978 ), is the resolution of the apparent con- 
flict between the luminosity ratio and the mass ratio. Our pho- 
tometric solution indeed yields a luminosity ratio close to 3.1. 
This is well below the value (by a factor of 100 or more!) ex- 
pected from the dynamical mass ratio of 3.3, if we were dealing 
with single stars or components of an unevolved detached binary 
system! We note also that the secondary's absorption lines yield 
an optical brightness ratio of 1.4 + 0.6 (Rauw et al. 1999) which 
is in marginal agreement with the value of 2.9 inferred from the 
photometric solution. 

The key to this puzzle probably lies in the fact that the two 
stars are in a contact configuration. In their pioneering work to 
explain the unusual mass-luminosity relations of W UMa-type 
contact binaries, Osaki (119651 ) and Lucy ( 119681 ) showed that 
the fact that energy is radiated according to a gravity darken- 
ing law implies that in a common envelope system ^ ~ it^ 
which agrees with our results for CygOB2 #5. These authors 
also showed that both stars in a contact binary system should 
have almost equal surface brightnesses, which is indeed what 
we observe in CygOB2 #5. Whilst Osaki ( 119651 ) adopted the 
conventional von Zeipel gravity darkening law that applies to 
stars with a radiative envelope, Lucy ( 1968 ) argued that the von 
Zeipel law is inappropriate for W UMa contact binaries that are 
typically of spectral type K - F and should thus have a convec- 
tive common envelope. Lucy (119681 ) also questionned the possi- 
bility to apply Osaki's (1965) conclusions to early-type contact 
binaries (such as the one we are dealing with in this work), be- 

' Concerning the latter, we note that INTEGRAL observations (De 
Becker et aI. I2007t failed to detect a hard X-ray, low y-ray counterpart 
to the EGRET source. 
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cause the common envelope would not be in hydrostatic equilib- 
rium. However, Tassoul (119921 ) showed that early-type contact 
binaries can have a barotropic outer common envelope hence 
validating the use of von Zeipel's gravity darkening law. Whilst 
Cyg OB2 #5 is certainly far hotter and more massive than classi- 
cal W UMa-type systems, it seems likely also in our case that the 
spectra of both components originate in a common photosphere 
and this explanation then accounts for many observed proper- 
ties of the system. However, there is one caveat; our analysis has 
revealed a secondary's star surface that is apparently not in equi- 
librium since we found a bright region facing the primary star 
Therefore, more sophisticated theoretical models are needed to 
understand the structure of a common envelope during mass ex- 
change in early-type binaries. 

In previous papers, it was argued that the secondary of 
Cyg OB2 #5 might be a Wolf-Rayet or Ofpe/WNL transition ob- 
ject with a spectrum altered by a hydrogen cloaking due to mass 
transfer from the primary (Bohannan & Conti 1976; Vreux il985i 
Rauw et al. 119991 ). In the light of the results obtained in this pa- 
per, we emphasize that there is no real need for such a scenario 
anymore. The contact configuration indeed probably accounts 
for most of the odd features of this star, which could otherwise 
be a rather 'normal' Bl-2 star 

4.4. The distance of Cyg 0B2 

Our analysis of the light curve of Cyg OB2 #5 yields a distance 
of 900 - 950 pc, roughly half of the previous distance estimates 
of Cyg OB2. While this result needs to be confirmed, it is inter- 
esting to note that this revision of the distance would have im- 
portant implications for our general understanding of the stars in 
Cyg OB2. First of all, the luminosities (and luminosity classes) 
of the O stars would have to be revised. A reduction in the lumi- 
nosities and thus also the mass-loss rates of the stars would then 
have a potential impact on our understanding of the radio light- 
curve of the non-thermal radio emitting O-stars such as Cyg OB2 
#8a and #9 (De Becker et al. l2006l Naze et al. |2008] |. In these 
binary systems with relatively short orbital periods, the wind in- 
teraction zone is rather deeply embedded in the winds of the 
stars. With current estimates of the mass-loss rates, the bulk of 
the synchrotron radio emission produced in this wind collision 
region would be absorbed by the winds themselves. A downward 
revision of the distance of these stars and their wind densities 
would imply a lower opacity in the radio domain and could help 
to solve this issue. 

5. Conclusion 

Using optical narrow-band photometry and spectroscopy, we 
have shown that Cyg 0B2 #5 is in a contact configuration, which 
has allowed us to solve one of the oldest mysteries associated 
with this system, i.e. the apparent discrepancy between its mass 
ratio and its luminosity ratio. The binary appears to be in a rather 
'normal' stage of (likely case-A) mass transfer in binary evolu- 
tion. 

However, whilst we have made progress in understanding the 
system, there remain a number of major issues. The most impor- 
tant ones are: 

• New theoretical models need to be designed to simulate ac- 
curately the structure of a common envelope in an early-type, 
contact binary system. 

• The origin of the period variations needs to be investigated 
by long-term photometric monitoring of the system. 



• The origin of the X-ray emission in this system remains un- 
certain as is the cause of the variations in the radio emission 
(Kennedy et al. 120091 ). The most likely explanations involve 
the existence of a fourth component in the system. However, 
to date no direct evidence that such a component exists has 
been obtained and again long term monitoring of the system 
is required. For instance, one needs to acquire optical spec- 
troscopy to determine orbital solutions that sample the 6.7 yr 
cycle. 

• And finally, the distance of Cyg 0B2 needs to be checked 
by studying more eclipsing binary systems in this important 
cluster 
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